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ABSTRACT 

Important clues to the chemical and dynamical history of elliptical galaxies are 
encoded in the abundances of heavy elements in the X-ray emitting plasma. We 
derive the hot ISM abundance pattern in inner (0 — 2.3i? e ) and outer (2.3 — 4.6i? e ) 
regions of NGC 4472 from analysis of Suzaku spectra, supported by analysis of co- 
spatial XMM-Newton spectra. The low background and relatively sharp spectral 
resolution of the Suzaku XIS detectors, combined with the high luminosity and 
temperature in NGC 4472, enable us to derive a particularly extensive abundance 
pattern that encompasses O, Ne, Mg, Al, Si, S, Ar, Ca, Fe, and Ni in both 
regions. We apply simple chemical evolution models to these data, and conclude 
that the abundances are best explained by a combination of a-element enhanced 
stellar mass loss and direct injection of Type la supernova (SNIa) ejecta. We 
thus confirm the inference, based on optical data, that the stars in elliptical 
galaxies have supersolar [a/Fe] ratios, but find that that the present-day SNIa 
rate is ~ 4 — 6 times lower than the standard value. We find SNIa yield sets that 
reproduce Ca and Ar, or Ni, but not all three simultaneously. The low abundance 
of O relative to Ne and Mg implies that standard core collapse nucleosynthesis 
models overproduce O by ~ 2. 
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Subject headings: galaxies: abundances, galaxies: elliptical and lenticular, galax- 
ies: individual (NGC 4472), galaxies: ISM 

1. Introduction 
1.1. Context 

Galaxy formation encompasses the processes of star formation and largescale dynamics 
- where dynamics is defined in its broadest sense to include the assembly of dark matter and 
baryons, as well as the exchange of mass and energy among the various galactic components 
and with the external environment. Star formation ultimately leads to the nuclear production 
of metals, while dynamical processes determine their destination. Therefore the abundance, 
abundance pattern, and location of metals in different galaxies in the local universe provide 
fundamental insights into the construction and development of galaxies.. 

Chemical evolution modeling quantifies this connection. Heavy elements are both stat- 
ically and explosively synthesized in evolving stars following episodes of star formation with 
abundances that reflects the total mass of stars formed and the initial mass function (IMF), 
as well as the production rate of Type la supernova (SNIa) progenitor binary stars. These 
metals are subsequently returned to the interstellar medium (ISM), and the enriched ISM 
may regenerate stars or escape into intergalactic space, depending on the level of supernova 
energy injection that accompanies the metal enrichment - and how efficiently that energy is 
channeled into outflow. These processes must depend on initial conditions and environment 
in such a way as to produce the diverse and evolving universe of galaxies that we observe 
from the Local Group out to the distant universe at redshifts of 6 and greater. One can 
illuminate the history in the ensemble of galaxies in a cluster by measuring abundances in 
the intracluster medium (ICM; see, e.g., Loewenstein 2006), and in individual galaxies by 
measuring stellar and interstellar abundances. 



It is evident that in giant elliptical galaxies, predominantly composed of stellar popula- 
tions passively evolving since the early universe, the star-gas cycle was greatly accelerated and 



concentrated in time relative to the Milky Way and other late-type galaxies ( 


PaDovic 


i et al. 


2006; 


van Dokkum & van der Mare! 


2007; 


Jimenez et al. 


2007; 


Perez- Gonzalez et al. 


2008h. 



Moreover, the prodigious metal content in the intergalactic medium of galaxy clusters domi- 
nated by ellipticals implies that strong galactic winds were driven in these galaxies. However, 
details such as the precise initial epoch and duration of star formation, the form of the IMF, 
the relative roles of SNIa and core collapse (SNII) supernovae in enriching and expelling ISM 
- and how all these depend on mass and environment - remain unclear. 
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Abundances in the stellar component and their correlation with structural parameters, 
provid e key diagnostics of elliptical galaxy ages, star formation and accretion histories, and 



IMFs (iThomas et al.l20 05: Pip ing et al 



2009 



Clemens et al.ll2009t IColavitti. Piping. & Matteucci 



20091 ; iTrager & Somervillell2009l ; iRecchi et al.ll2009l ) . Since stars dgminate the barygn cgntent 
gf ellipticals, the level gf stellar a-element enrichment is determined by the past integrated 
frequency gf SNII, and hence the number gf massive stars formed. Fgr a given total mass in 
stars, this then constrains the IMF. The stellar a/Fe ratio is then determined by the relative 
contributions of SNIa and SNII. Careful modeling of multiband photometry and very high 
signal-to-noise optical spectroscopy can, in principle, determine the abundances as well as 
the ages of ellipticals. However, since interpretation of colors and absorption line indices 
suffer from a degeneracy in their dependence on the S FH and overa ll metallicity - even for 
an assumed IMF and distribution of abundance ratios (jHowelll 120051 ) - optical studies of the 
composite stellar population cannot clearly and unambiguously measure these fundamental 
quantities. The additional co nsideration of Balm er emission line indices can assist in break- 
ing some of the degeneracies ( ITrager et al.l 120001 ) but, e.g., does not easily distinguish a very 
old stellar po pulation with a "frost i ng" of recent star formatio n episodes from a somewhat 
younger one ( jSerra fc Tragerl 120071 ; ITrager fc Somervilld 120091 ). Since, optical line indices 
can only be interpreted within the context of a conjectured SFH and IMF, they are more 
robust in providing consistency tests of combinations of ag es and abund ance ratios (i.e., 
Mg/Fe) than they are in determining absolute abundances ( jHowelll 120051). Moreove r, the 
effects of variations in individual elements is not easily disentangled (ILee et al.l 120091 ). and 
optical absorption features are only measured out to 1-2 (projected) effective radii (R e ) - 
the abundances in, typically, half of the stellar mass are unknown. 

Fortunately, mass lost by post-main-sequence stars in elliptical galaxies is heated to 
millions of degrees K as it is incorporated into the interstellar medium (IS M), where it is 



amen able to relatively robust and straightforward X-ray spectroscopic analysis ( iMathews fc Brighenti 



20031 ). Because the gas, with the possible exception of a few high oscillator strength lines, is 
optically thin the absolute abundances of a wide range of elements with prominent emission 
features that can include C, N, O, Ne, Mg, Al, Si, S, Ar, Ca, Fe, and Ni may be directly 
derived given sufficient spectral resolution, sensitivity, and bandpass. Because this approach 
relies on strong emission lines, X-ray derived abundances may be deter mined out to very 
large radii - often extending to the edge of the optical galaxy and beyond ( Matsushita et al. 
1998h . 



While hot ISM abundances may be more directly derived than stellar abundances, their 
interpretation requires careful deconstruction within the context of physical gasdynamical 
and chemical evolutionary models. Unlike the case of the ICM that dominates the baryon 
content in clusters, most of the metals produced by the evolving stellar population are 
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expelled from galaxies or locked up in stars - not accumulated in a reservoir of hot gas. 
Enrichment timescales for elements synthesized by SNIa, SNII, and intermediate mass stars 
are distinct, and coupled in a complex way to star formation, mass return, and outflow 
timescales. Nevertheless, one can employ reasonable assumptions, approximations, and sim- 
plifications to construct simple models that track the evolution of global abundances of these 
elements. Comparison of these models with observations constrain important features of the 
nature of the galaxy stellar population, as well as the rate and ultimate disposition of metals 
ejected as this population evolves. As such, they may serve as guides to subsequent more 
complex modeling. 

In this paper, we introduce such models and apply them to the hot ISM abundances in 
NGC 4472 that we infer from analysis of Suzaku and XMM-Newton spectra. Because of its 
brightness and temperature structure, a particularly well-determined and wide-ranging X-ray 
abundance pattern is measurable in NGC 4472 - one that is made more accessible utilizing 
the low internal background and sharp energy resolution of the Suzaku XIS detectors. 



1.2. Brief Survey of Previous X-ray Results 



High quality X-ray spectra for a large sample of elliptical galaxies was first obtained 
with the Advanced Satellite for Cosmology and Astrophysics (ASCA). Fe abundances in 
most of the X-ray brightest ellipticals were found to be within a factor of two of solar , 



with no strong evidence of non-solar r atios of a-elements to Fe (IBuote fe Fabian! 11998 
Matsushita. Ohashi. fe Makishimal 120011 ) . X-ray spectra extracted from the Chandra and 



XMM-Newton CCD detectors take advantage of improved sensitivity and spectral resolu- 
tion, and far superior angular resolution, to more cleanly derive abundance patterns and 
gradients (Humphrey & Buote 2006, and references therein). Chandra results indicate that 
ISM Fe abundances are roughly solar an d decline slowly with radius as far out as tens of kpc 
( [Humphrey fe Buotdl2006l ; lAtheyl 120071 ). In addition, approximately solar Mg/Fe and Si/Fe 
ratios and subsolar O/Fe ratios are inferred. 

Abundances derived from the XMM-Newton reflection grating spectrometer (RGS) ben- 
efit from cleaner separation of individual features, so that effects of resonance scattering on 
optically thick lines, and multiphase gas on temperature-sensitive lines and line ratios, may 
be directly addresse d. Analysis of X MM-Newton reflection grating spectrometer (RGS) spec- 
tra of NGC 4636 by lXu et al.l ( 120021 ) yields subsolar Mg/Fe and O/Fe ratios, a roughly solar 
Ne/Fe ratio, and a supersolar N/Fe ratio. The difficulty in analyzing grating spectra of 
extended sources, and limitations in sensitivity and bandpass, restrict the applicability of 
RGS spectroscopy to X-ray luminous galaxies with the highest central X-ray surface bright- 
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nesses. I Ji et al.l ( 120091 ) analyze Chandra and XMM-Newton EPIC and RGS data in 10 bright 
systems. They confirm the above abundance pattern with respect to O/Mg/Si/Fe, and ex- 
tend to Ne and S that show no strong deviations from solar ratios, and Ni that is generally 
supersolar. 

The lower internal background and sharper energy resolution of the Suzaku XIS CCDs 
enable the derivation of abundance patterns to larger radii, and allow for a more accurate 
measurement of possible features at high energy originating from S, Ar, Ca, and Ni. The 
Suzak u low energy sensitivity makes it s uitable for measurin g O as well. Results on NGC 



20091), NGC 507 f lSato et al. 



1399 (Matsushita et al.ll2007 ). NGC 720 (ITawara et al.l 120081). NGC 5044 (IKomivama et al. 



20091 ). and NGC 4636 (IHavashi et all 120091 ) confirm the solar 



Si/Fe and subsolar O/Fe ratios, demonstrate that these persist to large radii, and extend this 
behavior to S/Fe. It is the case that this is also broadly true of Mg/Fe; however, Mg/Fe is 
somewhat lower in NGC 720 and NGC 1399. Ne/Fe seemingly shows more variation, coming 
in at solar in NGC 4636, subsolar in NGC 720 and supersolar in NGC 507 and NGC 5044. 
Ne does not seem to trace O as might be expected. 

Optical spectra imply that the stellar ratio of a elements to Fe, expressed logarithmically 
with respect to solar as [a/Fe], is such that [a/Fe] ~ 0.3 — 0.5 - although there is little 
known about this ratio well beyond R e or about the relative abundances of different a 
elements. If stellar mass loss dominates ISM enrichment, they should display the same a- 
element overabundance. Based on the the above summary, this is clearly not the case. Direct 
injection of SNIa will skew that abundance pattern towards lower [a/Fe], but should do so 
in a way that does not effect relative abundances of elements with low SNIa nucleosynthetic 
yields such as O, Ne, and Mg. The ratios, with respect to Fe, of Si, S, Ar, and Ca may be 
less affected, depending on the level of SNIa enrichment since these elements are produced 
with moderate efficiency by SNIa. The solar ratio seen for most of the a elements in the 
hot ISM indicates a significant role for SNIa enrichment; but, the fact that O is particularly 
underabundant while Mg and Ne are as abundant as Si and S does not comfortably fit into 
this scheme. We examine these issues more quantitatively below with respect to the ISM 
abundances in NGC 4472, where measurements of Ca and Ar in the ISM add a new twist to 
attempts to understand the enrichment of this galaxy. 
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2. Construction of NGC 4472 Spectra and Associated Files 



2.1. Suzaku Spectral Extraction 



NGC 4472 was observed with Suzaku (IMitsuda et all 120071 ) between 2006-12-03 and 
2006-12-06 (OBSID=801064010) with an on-source exposure time of 121 ksec - 94 (27) 
ksec in 3 x 3 (5 x 5) CCD editing mode. At the time of observation, th ree co-aligned, 



17.8' x 17.8' field-of-view X-ray Imaging Spectrometer (XIS) CCD cameras (IKoyama et al. 



20071 ) - two front-illuminated (FI: XIS0 and XIS3) and one back-illuminated (BIXIS1) - 
were operational , each XIS in the focal plane of an X-ray Telescope (XRT) with a 2' half- 
power diameter (ISerlemitsos et al.l 120071 ) . 



Observations were conducted utilizing the space-row charge injection (SCI) technique 
that reverses the degradat ion in energy resolution caused by accumulated radiation dam- 
age (INakajima et al.ll2008l ). We initiate our data reduction with the unfiltered event files. 
These data underwent Version 2.0.6.13 pipeline processing on 2007-08-18 that enables one 
to properly account f or the effect of SCI on instrument characteristics and performance 
( jUchiyama et al.ll2008f ). We reprocess the unfiltered event files by hand in order to apply up- 
dated calibration data and software. Our analysis generally follows the procedures outlined in 
"The Suzaku Data Reduction Guide," □ as implemented in HEAsoft version 6.5.10 We recal- 
culate PI values and grades, select event grades (0, 2, 3, 4, 6) that correspond to X-ray photon 
events, filter on pixel status (eliminating bad charge transfer efficiency columns, and rows 
invalidated by the charge injection process) and select good time intervals (GTI) based on 
pointing, data and telemetry rates, SAA proximity ("SAA_HXD = 0, T_SAA_HXD > 436"), 
and proximity to the earth's limb and illuminated Earth ("ELV > 5, DYE_ELV > 20"). 
In addition, telemetry-saturated frames and calibration source photons are screened out; 
and, hot and flickering pixels are removed. Finally, we accept only GTI where the revised 
geomagnetic cut-o ff rigidity COR2> 4, thus eliminating intervals with the highest particle 
background level ( jTawa et al.ll2008l ) without compromising overall statistical accuracy (ex- 
periments with the stricter criteria COR2 > 6 yields identical results). 5x5 event files are 
converted to 3 x 3 mode format, and merged with the 3x3 event files. 

Spectra are extracted from — 4' (0 — 2.3R e ) circular (inner) and 4 — 8' (2.3 — 4.6R e ) 
annular (outer) regions, centered on the NGC 4472 optical nucleus (a = 12 29 46.7, 5 = 
+08 00 02) that very closely corresponds to the X-ray peak in the Suzaku image. Relatively 
large regions are chosen to minimize PSF effects, and maximize the accuracy and range 



1 http:/ /heasarc. gsfc.nasa.gov/docs/suzaku/analysis/abc/ 
2 http://heasarc. gsfc.nasa.gov/docs/software/lheasoft/ 
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of abundance measurements. Based on the simulations in ISato et al.l (120071 ). only a small 
fraction of counts in these regions originate in other annuli in the Abell 1060 galaxy cluster; 
and, this should be even more true for NGC 4472 based on its relative surface brightness. 

The spectral redistribution matrix files (rmf) are generated using xisrmfgen version 
2009-02-28. The rmf and spectral files are binned to 2048 channels. The effective area 
function files (arf) for the source spectr a are generated by the xissimarfgen version 2009- 



01-08 Monte Carlo ray-tracing program (jlshisaki et al.ll2007l ) with 400000 simulation photons 



per energy bin, and an input source fits image file generated from a /3-model fit to the 
background-subtracted X-ray surface brightness profile extracted from archival Chandra data 
(/3 = 0.44 and core radius 4", although the results proved insensitive in experiments with 
larger core radii ranging up to ~ 1' and/or = 0.48). Spectra from the FI chips, XIS0 
and XIS3, are co-added and a weighted XISO+3 response functions calculated from their 
respective rmf and arf files. Source spectra are grouped into bins with a minimum of 15 cts. 
Best-fit parameters and parameter uncertainties are derived using x 2 statistics; experiments 
with the Cash statistic yielded very similar results. 

The XIS background includes contributions from Non- X-ray (charged particle) Back- 
ground (NXB), Galactic X-ray Background (GXB), and (extragalactic) Cosmic X-ray Back- 
ground (CXB). Since NGC 4472 fills the Suzaku field of view, we estimate and subtract the 
NXB component and include an additional thermal component in our spectral fits to account 
for the GXB. The CXB has a similar shape, but a much smaller magnitude in the inner re- 
gion, than the intrinsic flux from low-mass X-ray binaries (LMXBs) and is subsumed into 
that spectral component. The NXB component is estimated from observations of the night 
earth taken in SCI mode within 150 days of the starting or ending dates of our observation 
using xisnxbgen version 2008-03-08. The NXB event list in that time interval undergoes 
the identical screening as the source data, is sorted by geomagnetic cut-off r igidity, and 



weigh ted according to the cut-off rigidity distribution in the source event file (ITawa et al. 



20081 ). The estimated NXB spectra include only those events collected in the regions on the 
detector from which the source spectra are extracted. Since the NXB has a spatial distri- 
bution distinct from that of the NGC 4472 X-ray emission, a separate arf file is generated 
using 2000000 simulation photons per energy bin from a uniform source of radius 20' and 
applied to the background in spectral fits. Combined XISO+3 NXB spectra and response 
functions are constructed as described above for the source spectra. 

The effective exposure time of the XISO+3 (XIS1) spectrum is 232 ks (117 ks). The 
total 0.3-12 keV band counts are 132032 (inner) and 70286 (outer) for XISO+3, and 107232 
(inner) and 72758 (outer) for XIS1. The fraction of the total count rate in the NXB is 1.4% 
(inner) and 8.3% (outer) for XISO+3 , and 1.8% (inner) and 10.3% (outer) for XIS1. Even 
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in the outer region, its contribution is small in the region where spectral features of interest 
lie (see Figure 3, below). 



2.2. XMM-Newton Spectral Extraction 

NGC 4472 was observed for ~90 ksec with XMM-Newton. After cleaning the dataset 
(the last 20 ksec of the full-field light curve showed significant flares), ~81 ksec (~72 ksec) 
of useful observation time remained for the EPIC MOS1 and MOS2 (pn) detectors. The 
data were reduced and extracted using SAS 6.5.0. The latest calibration products were used 
to determine the response files for the spectral analysis. The extracted spectra are binned 
so that each channel has a minimum of 20 counts, and we consider energy ranges of 0.4-7.0 
keV and 0.32-2.0 keV for the EPIC cameras and RGS detectors, respectively. To account 
for the fact that NGC 4472 is an extended source, we used the rgsxsrc model in Xspec 11, 
extracting the Chandra ACIS-S image to a create a surface brightness distribution model 
to convolve with the RGS response. The background for the central region was taken from 
the blank sky background files provided by the XMM-Newton Science Center. These were 
scaled to match the EPIC and RGS high energy continua. EPIC MOS1, MOS2, and pn 
spectra are extracted from the inner and outer regions described above. Subdividing the 
inner region (r < 200'') into multiple annular regions, and fitting the corresponding EPIC 
spectra, reveals no significant ~arcminute gradients in abundance ratios (Section 3.2) 



3. Spectral Analysis 

3.1. Suzaku Spectral Analysis 

The NGC 4472 Suzaku XISO+3 and XIS1 spectra are simultaneously fit in the 0.45-7.0 
keV band using Xspec 12.5 with a source model consisting of a two-temperature thermal 
plasma (vapec) model with a single set of heavy element abundances (C, N, O, Ne, Mg, 
Al, Si, S, Ar, Ca, Fe, Ni), plus a 7 keV thermal bremsstrahlung component to account for 
the LMXBs (and the emission from the CXB that is relatively insignificant in the inner 
region, and - for the inner region - any activity associated with the dormant supermassive 
black hole; see Loewenste in et al. 2001). The abs orption is fixed at the Galactic column 



density of 1.65 10 20 cm 2 (IDickey fe Lockmanlll990l ). When free to vary, the bremsstrahlung 



temperature and column density are found to be consistent with these values; and, the vapec 
parameters essentially unchanged. This is also the case if a power-law model is adopted to 
characterize the LMXBs. The inner region bremsstrahlung flux is consistent with previous 
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estimates of the LMXB contribution (lAtheyl 120071 ). and the outer region bremsstrahlung 



flux consistent with the expected sum of LMXB and CXB components. The Si K-edge 
region (1.82-1.841 keV) was excluded from the spectra, although this had negligible effect 
on the final fits aside from a small reduction in \ 2 - A const ant XIS1/XIS0+3 multipli cative 



factor is included and generally ranges from 1.01-1.05. The iGrevesse fe Sauvall (119981 ) solar 
abundance standard is adopted. Since no elements with features of significant strength are 
tied to solar ratios, the resulting abundances may be re-scaled to other solar standards. Of 
course, we use the same solar standard in the models that we compare with the results of 
spectral fitting. C emission lines are outside of the Suzaku bandpass, so the C abundance 
is fixed at solar. The inner and outer regions are separately fitted. The outer spectrum 
requires an additional cooler component, presumably from GXB emission that we model by 
a single temperature thermal plasma model with all abundances except N and O fixed at solar 
(allowing other elements to vary, or adding a second temperature component did not improve 
fits or affect hot ISM parameters of interest). Because the GXB is more prominent there, the 
outer spectra are used to determine the GXB parameters (except for the normalization) that 
are then fixed in analyzing the inner spectra. The relative fluxes of the GXB component 
in the inner and outer regions are consistent with the ratio of their projected areas, and 
the overall GXB surface brightness, temperature (0.15 ± 0.04 keV), and oxygen abundance 
(0.15, 0.11-0.29 at 90% confidence) are consistent with what was found in th e source-free 



Ursa Minor Dwarf Spheroidal Field (ILoewenstein. Kusenko. fe Biermannll2009l ) - there is no 



evidence that any of the "GXB" constitutes an intrinsic NGC 4472 cool component. Fits to 
the inner spectra require two ISM temperature components. The temperatures we find in 
the inner region are 0.798 ±0.004 and 1.42 ±0.05 keV. The outer region may be fitted with a 
single-temperature (1.2 keV) model (reduced \ 2 — 1-05), but is better fitted (A% 2 = 96) with 
a two-temperature (0.93 ± 0.08 and 1.55 ± 0.10 keV) model. These are consistent with the 
temperature structure we derive from spectra extracted over a finer series of annuli (see next 
section). The best-fit N abundances are with upper limits of 0.75 (0.6) solar in spectral 
fits to the inner (outer) spectra, and are fixed in the final fits. The addition of a third ISM 
temperature component did not significantly affect the best-fit abundances. 

Spectra (with residuals to the best-fit model) and best-fit models, are displayed in detail 
in Figures 1-4. The best fit abundances and 90% (Ax 2 = 2.7) confidence limits are shown in 
Table 1. The best-fit abundance ratios, with respect to Fe, for the inner and outer Suzaku 
spectra are compared in Figure 5. Note that O, Mg, Si, S, Ar, Ca, and Al (K-) line features 
are well-resolved (though the latter is weak), but that Ne (K-) and Ni (L-) line features are 
not, with respect to the Fe L-line complex. 
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Fig. 1. — top (a): XISO+3 inner region spectrum and best-fit model (black data points 
and histogram). Also separately shown are the GXB thermal plasma (orange), thermal 
bremsstrahlung LMXB (blue), and hotter (green) and cooler (red) hot ISM thermal plasma 
components, bottom (b): Same for XIS1. 
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Fig. 2. — Same as Figure 1, zoomed to the Si/S/Ar/Ca spectral region. Prominent features 
from left to right predominantly originate in the following ions: Si XIII (Ka), Si XIV {Ka), 
Si XII (K/3), S XV (Ka), S XVI (Ka), S XV (K/3), At XVII (Ka), Ca XIX (Ka). 
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Fig. 3. — Same as Figure 1 for outer region spectra, with NXB (purple) contribution and 
sum of non-background components (orange) now included. 
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Fig. 4. — Si/S/Ar/Ca spectral region as in Figure 2 for outer spatial region, including NXB 
(purple) contribution. 
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Fig. 5. — Inner (crosses), and outer (open triangles, offset to right) Suzaku abundance 
patterns. Ca is shifted to the right. 
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3.2. XMM-Newton Spectral Analysis and Comparison to with Suzaku 

XMM-Newton spectra are analyzed using the same models that we applied to the Suzaku 
analysis. Fits to the XMM-Newton EPIC M0S1+2 and pn data in the same extraction 
regions yielded best fits that are generally consistent with the above, though with some 
discrepancies (Table 1, Figure 6). Although blank sky background is subtracted from these 
spectra, an additional small GXB component is included. The best-fit GXB normalization is 
3 — 10 x lower than in the (non-GXB-subtracted) Suzaku spectra. The effect of this addition 
is to improve the spectral fits (A% 2 = 53 for the MOS, 16 for the pn) and increase the 
parameter uncertainties. As with the Suzaku spectra, a two-temperature plasma is required 
for the hot ISM. 

There are substantial differences in the derived best-fit abundances of Ar, Ca, and Ni 
derived from EPIC and Suzaku spectra (Figure 6a); but, the errors are large (particularly for 
the former) and the 90% confidence ranges overlap. The most serious discrepancies are the 
low outer-region Mg abundance and inner-region O abundance derived from EPIC spectra. 

The Suzaku abundance pattern we derive in the inner region is in very good accord 
with that found by Ji et al. (2009) from analysis of EPIC and RGS data, and from our 
own RGS analysis (Figure 6a, Table 1). Although the latter samples a smaller region, the 
best-fit temperatures and abundances are consistent with our analysis of the inner region 
Suzaku data except for N (best RGS fit: 3.7 ± 1.7 at 90% confidence) that is sensitive to the 
treatment of the background, and (marginally) Mg. 

The XMM-Newton EPIC angular resolution allows us to determine the temperature and 
abundances in 9 annular rings out to 300". In our analysis we find that a two-temperature 
vapec plasma model is needed for most of the radial bins. This requirement is driven by 
residuals to one-temperature fit s in the CCD spectra near the energy of the Fe XVII line. This 



is in contrast to the analysis of iNagino fe Matsushita ( 120091 ). who adopt a one-temperature 



model in deriving a deprojected temperature profile. In our analysis we have used fairly wide 
annuli (30" minimum) and adjusted their widths to include at least 20000 counts in order 
to allow us to accurately determine the abundances in each ring. This requirement results 
in extraction regions that span a larger fraction of the temperature gradient so that two- 
temperature models are needed to fit the data. Based on a comparison with one temperature 
models in the innermost annulus, the two temperature model is required at the 99.99% level. 
All other annuli require a two temperature fit with at least a 90% confidence, except for the 
outermost ring where the two temperature fit is only significant at the 80% level. We adopt 
the two-temperature model for this annulus as well to facilitate comparison with the other 
annuli. Our temperature profile is shown in Figure 7, and shows a positive gradient. The 
cooler component temperature is about 60% that of the hotter at all radii. The ratio of hot 
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to cool components increases with radius, so that in the outer annuli the hotter component 
dominates the emission. 

Utilizing these same annular extraction regions, we search for abundance gradients on a 
finer scale than that of the Suzaku data for the best-determined abundances: O, Ne, Mg, Si, 
S, and Fe. We determine the abundance gradient for each element with respect to Fe; and, 
in general the profile is consistent with a constant. The only ratio that shows a gradient is 
Mg-to-Fe, which is flat in the inner 150 - 200" and declines to about half the central value 
at ~ 300". 

In the following sections, we adopt the abundances that we infer from fitting the Suzaku 
spectra for comparison with simple chemical evolution models. We consider these more 
robust, for the crude spatial binning that we employ, because of the lower background and 
sharper energy resolution of Suzaku. This is illustrated in Figure 8 that compares MOS1, pn, 
and XIS0 inner region spectra, and Figure 9 that compares pn and XIS0 outer region spectra. 
The Suzaku spectral features, including the 1.0 keV Ne X bump embedded in the Fe-L line 
complex (note, also, the excellent agreement of the Suzaku inner spectrum Ne abundance 
with that derived from the RGS spectrum where this feature is resolved), the OVII line at 
0.65 keV, and the Mg XI and Mg XII lines at 1.35 and 1.47 keV are significantly sharper 
than their EPIC counterparts (and in some instances sharper than the resolution in the 
latest Suzaku calibration). The reduced high energy background is particularly significant 
in the Si/S/Ar/Ca spectral regime for the outer region (Figure 8). 

We note the formal measurement of a large Al abundance in the inner region (with 
XMM-Newton and Suzaku) and the outer regions {Suzaku only - there was residual non- 
X-ray Al fluorescence emission in the outer-region MOS spectra). This would be the first 
detection of an odd-atomic-number element in the hot ISM of an elliptical galaxy. An Al 
abundance tied at solar results in a poorer fit over the 1.5-1.8 keV region (A% 2 = 29 for the 
inner Suzaku spectrum). However, the inferred 1.6 keV Al XII and 1.7 keV Al XIII features 
are weak, and possibly sensitive to uncertainties in the characterization of the continuum in 
their vicinity. We do not include Al abundances in the model constraints discussed in the 
subsequent sections, but do sometimes display them in comparisons of models with data. 



4. Interpreting the NGC 4472 Abundance Pattern 



Our method of interpretation, described in more technical detail in Appendix A, may be 
summarized as follows. We assume that the ISM abundances are dete rmined following the 
cessation of the main era of star formation < 9 Gyr ago for NGC 4472 ( [Humphrey fe Buote 
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Table 1. Best Fit ISM Abundances 
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Note. — Reduced \ 2 an d best fit abundances relative to the Grevesse & Sauval (1998) solar standard, as well as 90% confidence 
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Fig. 6. — left (a): Inner region abundance pattern (normalized to Fe) derived from Suzaku 
XISO +3 and XIS1 (crosses), XMM-Newton EPIC from our analysis (open squares) and 
from lJi et all ( 120091 ) (open triangles), and XMM-Newton RGS (open polygons), right (b): 
Outer region abundance pattern (normalized to Fe) derived from Suzaku XISO+3 and XIS1 
(crosses), and XMM-Newton EPIC (open squares). 
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Fig. 7. — The projected temperature profile for N4472 from the center out to 300". The 
crosses denote the low temperature component, and the filled circles the hotter component. 
The error bars show the 90% confidence range. The horizontal bars denotes the radial range 
over which the ISM spectral parameters were derived. 
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Fig. 8. — top (a): MOS1 (red), pn (green), and XISO (black) inner region spectra in the 
0.6-1.7 keV bandpass, bottom (b): Same for 1.75-4 keV bandpass. 
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Fig. 9. — top (a): pn (red), and XISO (black) outer region spectra in the 0.6-1.7 keV 
bandpass, bottom (b): Same for 1.75-4 keV bandpass. 
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20061 ; iGallagher et all 120081 ). at which time the stellar abundances are established. The 



internal sources of mass and metals are stellar mass loss - calculated using remnant masses 
and main sequence turnoff times from the literature - and SNIa explosions. Effects of outflow 
(including the possibility of a mismatch between the metallicity of the gas in the galaxy and 
that flowing out of the galaxy), initial conditions, and different time dependencies for the 
mass return and SNIa rates are considered. For a set of simple one-zone models ( "closed- 
box" , and "tracking" solutions where mass return and SNIa have the same time-dependence 
- steady state solutions are a subset of the latter), we solve the equations expressing the 
conservation of overall ISM mass and mass of the ith element for the present-day ISM 
abundances. For these simple models, solutions generally are of the form 

/iSM = r /stars + S( ?Z/SNIa) (1) 

where /j SM and / s \ ars are the ith element ISM and stellar mass fractions, respectively; 
r is a stellar mass loss dilution (r < 1) or enhancement (r > 1) factor, and sq is a 
term that determines the importance of SNIa enrichment relative to that from stellar mass 
loss. As explained in Sections Al-3, q is the ratio of present-day SNIa to mass return 
rate, q = 0.006676 ) sNia/6 ) MR Mq 1 , where ^mr is the specific mass loss rate in units of 
2.4 10~ n M Q yr _1 L~j^ & (w e henceforth assume # mr = 1), and ^sNia is the specific SNIa 
rate in units of 0.16 SNU (jCappellaro et al.lll997l ). The dimensionless r and s parameters 
may be element-dependent, but we do not consider that generalization here. The input 
parameters are (1) the stellar [ct/Fe] ratio that we use to determine the stellar abundance 
pattern usi ng (IMF-weighted, so la r metallicity progenitor ) SNII+HN (hypernova) and SNIa 
yields from Nomoto et al. ( 1997 ); Kobayashi et al. ( 20061 ) . (2) the ratios of the SNIa yields 
for each element relative to that of Fe, (3) the ratio of ISM Fe abundance to stellar Fe 
abundance, and (4) the dilution parameter r. The output is the ISM abundance pattern 
relative to Fe. Specifying the ISM and stellar Fe abundances and SNIa yields (not just 
the ratios), determines sq and the predicted values of the ISM abundances. The r and sq 
parameters can then be interpreted in terms of relative SNIa frequency, outflow parame- 
ters, and initial condition in the context of different models (Section A3). We generally 
consider [a/Fe] sta rs = 0.25,0.15,0,-0.2 (where the sum used to compute a includes only 
those elements measured in the ISM), an ISM/stellar Fe abundance ratio 1.65 based on a 
gas-mass weighted average of our inner a nd ou ter region Suzaku abundances and the stellar 
Fe abundance from iHumphrey fe Buotd (120061 ). and r = 1, 2/3, 1/3. 

A value of r < 1 may be explained in several ways (Sections A.1-A.3). In the steady- 
state models we consider, r = 1/3 (2/3) corresponds to an outflow bias of 2.0 (0.5) - i.e, that 
gas is outflowing with 1.5 x (3.0 x) the mean metallicity inside the galaxy; r = 1 corresponds 
to no such bias. For the closed box models with the extreme initial condition of initial ISM 
metallicity, r = 1/3 (2/3) corresponds to having two-thirds (one-third) of the present-day 
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ISM mass at 84% (91%) of the Hubble time. The dilution is attributed to this preexisting 
low-metallicity gas; r = 1 corresponds to (the reasonable) initial conditions where the ISM 
and stellar abundances are equal. For tracking solutions, where winds are driven by SNIa 
and the mass return and SNIa rates have identical time-dependencies, we can also consider 
the extreme case where the ISM metallicity is at the initial time of 0.7 of the Hubble time. 
Here, r = 1/3 (2/3) corresponds to having 1.8 x (35 x) the current ISM mass. The dilution 
is due to the large initial reservoir of low-metallicity gas that is driven out of the galaxy 
while being replenished by SNIa-enriched stellar mass return; r = 1 again corresponds to 
(the reasonable) initial conditions where the ISM and stellar abundances are equal. 



We consider two sets of SNIa yields from iNomoto et al.l ( 119971 ) - the standard W7 de- 



flagration model, and the WDD1 delayed detonation model that represents the most radical 
departure from W7 in terms of nucleosynthetic yields. A comparison of the Suzaku-derived 
ISM abundance pattern inside the optical galaxy (the gas-mass weighted sum of our inner 
and outer abundances) and the implied stellar abundance pattern is shown in Figure 10. For 
stellar abundances derived using the W7 yields, the ISM abundances are similar to those in 
the stars assuming a solar or slightly supersolar [a/Fe} sta , rs - with the notable exceptions of 
ISM overabundances of Ar and Ca and an underabundance of O (all relative to Fe). For 
the WDD1 yields, the ISM again is similar to the [a/Fe] Btais ~ pattern, but now Ni is 
overabundant while O and S are underabundant in the ISM. Since direct SNIa enrichment is 
expected to distort the ISM pattern relative to that in the stars, these comparisons should 
be considered referential; however, some of these anomalies persist. 

To review how we compare these models with the observed abundances, sufficient SNIa 
enrichment (value of sq) is added to the diluted (since we consider r < 1) stellar Fe abun- 
dance to match the observed ISM Fe abundance. The SNIa rate corresponding to "sufficient 
enrichment" is model-dependent - but only weakly so for quiescent systems, where s ~ r, 
i.e. q pa (sq)/r or 6*sNia ~ 150(sg)/r. The magnitudes of the stellar mass loss and SNIa 
terms are now set and may be used, in conjunction with the stellar abundance pattern and 
SNIa yields, to determine the entire ISM pattern. Figure 11 compares model predictions of 
the ISM abundances with the observations for the case where W7 yields are adopted; #sNia 
corresponding to s = r is indicated^ The abundances of Ne, Mg, Si, S and Ni (as well as 
Al) are well-explained as a combination of undiluted (r = 1) stellar mass loss from stars 
with [a/Fe] s tars=0.25 and SNIa enrichment at a rate approximately one-sixth the standard 



3 We note here that our results are not sensitive to our adoption of the Grevesse fc Sauval ( 1998 ) solar 
abundance standard, since all quantities are scaled on this basis. Adoption of a different standard would 
shift observed and model abundances in lockstep, although the value of [a/-Fe] s tars corresponding to a given 
model would shift somewhat. 
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Fig. 10. — Comparison of abundances (including Al - located between Mg and S), relative 
to Fe, in the ISM (crosses with errorbars) and that expect ed in the stars for th e cases where 
W7 SNIa (a, left) and WDD1 (b, ri ght) SNIa yields (iNomoto et al.l [l997h are used in 
conjunction with SNII+HN yields from Kobayashi et al. ( 20061 ) to set the stellar abundance 



pattern for a given [a/Fe] Btaia . Shown are [a/Fe] sta rs=0.25 (solid line), 0.15 (dotted), 0.0 
(short-dashed) and -0.2 (long-dashed). 



-24- 



early-type galaxy value, #sNia = 0.17 (Figure 11a). However, Ar and Ca are underpredicted 
by factors of order 2 and O overpredicted by a similar factor. On the other hand the model 
with sufficient dilution to imply #sNia = 1 for a steady-state model (r = 1/3), the expected 
low ratios with respect to Fe of Ne, Mg, Si, S, Ar, Ca are greatly at odds with the obser- 
vations; and, the Ni/Fe ratio is overpredicted (Figure lib). In the limit where the stellar 
abundances are purely determined by SNII ([a/Fe] stars ~ 0.43), the ISM O-to-Fe ratio is 
reproduced but all other ratios with respect to Fe are not (Figure lib). 

The WDD1 models, where SNIa produce more Si, S, Ar, and Ca and less Fe and Ni, 
r = 1 models with supersolar [a/Fe\ stairs now adequately explain the ratios, with respect to 
Fe, of Ne, Mg, Ar, and Ca (as well as Al) but overpredict Si (slightly) and S as well as 
(Figure 12a). Ni is underpredicted by ~ 2, as it is for the r = 1/3 model. In the latter 
(Figure 12b), Si and O are matched by models with [a/ Fe] staxs corresponding to a pure SNII 
origin, but Ne and Mg are now underpredicted. 



5. Discussion 

In the most natural models - whether they resemble the closed box, steady state, or 
slightly more complex outflow models considered here - abundances (expressed as mass 
fractions) in the hot ISM of elliptical galaxies should be well-approximated by 

/bm(1) = fLs + 0.00667fl S NIa2/sNIa- (2) 

Figures 11 and 12 clearly show that the modest overabundance of interstellar Fe with respect 
to that in the stars, and the nearly solar abundance ratios of Ne, Mg, Si, and S with respect 
to Fe in NGC 4472, are consistent with this expectation provided that #sNia ~ 0.2 and 
[a/Fe] s tars ~ 0.2. We thus confirm the inference, based on optical data, that the stars in 
elliptical galaxies have supersolar [a/Fe] ratios - and infer that this extends globally. The 
relative Fe abundances may be accommodated in models with #sNia ~ 1 only by introducing 
mechanisms such as preferential outflow of metal rich gas, or assuming that the ISM has 
built up its Fe abundance relatively recently from substellar values either as the ISM mass 
itself has increased (for a closed box), or in concert with strong outflows. However, the ISM 
abundance pattern is poorly reproduced in such models for virtually all elements measured 
in the NGC 4472 ISM irrespective of the assumed [a/Fe] s t ars . 

While we consider only a single galaxy, large ISM-to-stellar Fe abundance ratios have 
not been measured in any elliptical galaxy (see references in Section 1.2). There is thus a 



conflict between the canonical elliptical galaxy SNIa rate (jCappellaro et al.l 119971 ) and the 



Fe abundance in the hot ISM of ellipticals. One avenue of reconciliation is to presume that 
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Fig. 11. — Comparison of abundances (including Al - located between Mg and S), relative to 
Fe, in the ISM (crosses with errorbars) and those predicted in models where W7 SNIa yields 
are adopted for r = 1 (left a:) and r = 1/3 (right b:). Shown are [a/Fe] stars =0.25 
(solid line), 0.15 (dotted), 0.0 (short-dashed) and -0.2 (long-dashed). Additionally, the 
[a/Fe] s tars = 0.45 (pure SNII stellar abundances) model is shown (dot-dashed line). 
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Same as Figure 11 for WDD1 SNIa yields. 



-26 - 



the effective rate of SNIa enrichment is much less than the actual rate because SNIa are not 
efficiently well-mixed into the ISM. This implies that there is large amount of missing Fe. 

It is plausible that there is substantial galaxy-to-galaxy variation in the current SNIa 
rates of elliptical galaxies. Galax ies that are radio loud or have blue colors have been 
shown to have enhanced SNIa rates (IMannucci et al.ll2005l ; iDella Valle et al.ll2005i ). Currently 
quiescent galaxies may be sampling the tail-end of a SNIa delay-time distribution with a peak 
at one or more timescales much less than the Hubble time, and the rates may be sensitive to 
the strength and lookback time of minor star formation episodes. If this is the case, galactic 
winds might be induced in those galaxies with the highest rates, rendering them gas poor 
and X-ray dim, and creating an observational bias tilted towards ellipticals with low SNIa 
rates. 

There are robust anomalies in the abundance pattern that cannot be explained using 
any of the standard yields (Figures 11 and 12), and this is the case even if we narrow our 
consideration to those elements with spectral features well-resolved by Suzaku. The ISM 
abundances of O, Ne, and Mg are overwhelmingly determined by SNII ejecta incorporated 
into, and then ejected, by evolved stars. Therefore, the NGC 4472 O underabundance 



universally rep licated in other elliptica l galaxy hot ( Humphrey fe Buotd 120061 ; I Ji et al. 



20091 ) and warm (lAthey fe Bregmanll2009l ) ISM - can only be explained by decreasing the O 
abundance in mass losing stars that, in turn, implies a reduction in (IMF-weighted average) 
SNII yield by about a factor of 2. This is not physically unreasonable considering the 
sensitivity of SNII production of O to the treatment of convection in the progenitor stars 
( iGibson. Loewenstein. fe Mushotzkyi 119971 ) . Alternatively, one could reduce the predicted 
O yields by a truncation of the IMF, since O yields increase steeply with progenitor mass; 
however, the mass dependence of Ne and Mg yields are nearly as steep (e.g., Loewenstein and 
Mushotzky 1996). These results for O are especially significant since the [O/H] abundance 
in HII regi ons is often used as a proxy for overall met allicity in studies of galaxies with star 
formation (lGarnettll2002l ; iKobulnicky fe Kewleyll2004l ). Interpretation of such data assuming 
solar ratios among a-elements and standard O yields should be considered with caution, as 
should calibration of SN-induced feedback in semi-analytic models that are based on observed 
O abundances. 

Another evident abundance pattern anomaly manifests as simultaneously high ISM 
abundances of Ar, Ca, and Ni. While application of W7 models can explain the high Ni/Fe 
ratio, it underpredicts Ar and Ca (Figure 11); the converse is true for models utilizing the 
WDD1 SNIa yields. For the r = 1, [a/Fe] stars =0.15 and 0.25 models with W7 SNIa yields, 
stellar mass loss accounts for > 80% of O, Ne, Mg, Si, S, Ar, and Ca, ~ 60% of Fe, and ~ 50% 
of Ni in the hot ISM with the remaining metals originating from direct enrichment via SNIa. 
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To "fix" the Ca and Ar discrepancy by increasing the stellar mass loss abundances through 
a boost in the SNII Ar and Ca yields is not ve ry appealing, since Ca is not overabun dant 



(with respect to Mg) in red-sequence galaxies ([Graves et al. 



2007; 



Smith et all 120091 ) - in 



fact, Ca traces Fe more closely than it does Mg ([Graves et al.l 120071 ). which argues for an 
important SNIa contribution to Ca. In Table 2 and Figure 13, we show the output for an 
r = 1 (6*sNia ~ 0.25), [a/Fe] st ars=0.15 model with the following ad hoc yield adjustments: a 
decrease in the average SNII O yield by a factor 2, and SNIa yields for Ca, Ar, Fe, and Ni 
of 0.05, 0.05, 0.5, and 0.07 M Q , respectively. One can see that, in such models, direct SNIa 
enrichment and enrichment via stellar mass return from material originating in SNIa and 
in SNII all make significant contributions to these elements, while the abundances of the 
lower-Z elements predominately originate in SNII products incorporated into stars. This is 
not meant as a definitive proposal, but to illustrate the sort of adjustments required. Should 
such abundance pattern anomalies prove prevalent, genuine constraints on supernova yields 
may be derived (Loewenstein, et al. in preparation). 



6. Summary and Conclusions 

We have derived the hot ISM abundance pattern in inner (0 — 2.3R e ) and outer (2.3 — 
4.6i? e ) regions of NGC 4472 from analysis of Suzaku spectra and used them to study the 
chemical evolution in this elliptical galaxy. The low Suzaku background and relatively sharp 
spectral resolution of the Suzaku XIS detectors enabled us to extend the range of accurately 
measured abundances beyond what is feasible with Chandra or XMM-Newton to S, Ar, 
Ca, and possibly Al - while we find general agreement with overlapping XMM-Newton 
results. The abundances of Ne, Mg, Si, S, Fe, and Ni may be explained by enrichment via a 
combination of a-element enhanced stellar mass loss and direct injection of SNIa with W7 



Table 2. Ad hoc Model Yields and Origins 
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Fig. 13. — Comparison of abundances (including Al - located between Mg and S), relative 
to Fe, in the ISM (crosses with errorbars) and those predicted in a model where ad hoc yields 
for some elements (see text, Table 2) are adopted for r = 1, [a;/Fe] stars =0.15 (6*sNia ~ 0.25). 
The latter is indicated by the solid line; stellar abundances are shown by the dotted line. 
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yields exploding at a current rate of 0.027 SNU; however, additional SNIa production of Ca 
and Ar is required. Ca and Ar are reproduced in models (with a SNIa rate of 0.035 SNU) 
where WDD1 yields are adopted, but at the price of disrupting the previous concordance for 
S and Ni. 

Models with the standard SNIa rate of 0.16 SNU badly overpredict Fe, unless the en- 
richment from stellar mass loss is diluted. We introduced models where this is a consequence 
of outflow biased in favor of high metallicity, or an inertial effect resulting from the preex- 
istence of a metal poor ISM. However, such models predict an underabundance of elements 
primarily synthesized in SNII (Ne and Mg for W7 and WDD1 SNIa yields, and Si, S, Ar, 
and Ca as well for W7) and offer a poorer match to the observed abundance pattern. 

We confirm previous measurements of a low O abundance; the low O/Mg and O/Ne 
ratios imply that the SNII yields in standard models overproduce O by ~ 2. 

Analysis of additional Suzaku and XMM-Newton spectra of elliptical galaxies will be 
used to investigate the universality of these results, and probe for diversity in the star for- 
mation and SNIa histories of giant elliptical galaxies. Application of more rigorous chemical 
evolution models will sharpen these conclusions, and should yield new information on the 
formation and evolution of these systems and indirect constraints on the physics of Type la 
super novae. 

A. Simple Models for Elliptical Galaxy Hot ISM Abundances 

A.l. Formalism 

In this paper we apply a simplified chemical evolution formalism to the NGC 4472 
observations. A more comprehensive treatment and application of the coupled chemical 
evolution of the gas and stars in elliptical galaxies is in progress. 

For a single-phase interstellar medium in a passively evolving galaxy - good approxima- 
tions for giant elliptical galaxies in general, and NGC 4472 in particular - the equations for 
the conservation of overall ISM mass, M ISM , and mass of the ith element, /i SM M ISM where 
/ism 1S ^e mass fraction, are as follows: 



and 



(Al) 
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d(fhuMisu) 
dt 



'-'SNIa u 



MR 



out J ISM" 



(A2) 



The only significant internal source of gas is stellar mass return (injected at the rate 
Mmr), while metals are provided by Type la supernovae (rate Gg NIa ) as well as stellar mass 
loss (rate G l MR ). M out accounts for gas flowing in or out of the optical galaxy (and is positive 
for outflow). The "bias" factor b l ont allows for an offset between the in- or outflowing gas 
metallicity and the mass-weighted-average inside the galaxy; henceforth we assume this is 
identical for all elements and suppress the index. For b out > (<)0 winds serve to reduce 
(increase) the overall galactic ISM metallicity, while for 6 out = winds remove metals but 
do not change the overall metallicity. For inflow, b out > (<)0 increases (reduces) the galactic 
ISM metallicity, 

These equations may be integrated from the end of the galaxy formation era t; n (at 
which time it is assumed that the galaxy is fully assembled and star formation has ceased) 
to the present day t now (13.7 Gyr). It is sensible to "start the clock" at t in since we know 
that the ISM mass inside the optical galaxy today is much less than the stellar mass return 
rate integrated over the galaxy lifetime. The star formation era must have been accom- 
panied by prodigious outflows in a multiphase I SM, a nd cannot be modeled using the for- 



malism developed here. Following ICiotti et al.l ( 



991), and ut ilizing main s e quenc e turnoff 



(11992( 1 and ISalaris et al.l (120091 ) . respec- 



masses and remnant masses from ISchaller et al 
tively, we derive Mmr ~ AL B {t/t now )~ a , where a fa 1.3, A fa 2.4 10 _11 #mr MqJt^ 1 L 



BQ 



and 6*mr = 1 for a simple stellar population with present-day blue luminosity Lb- In this 
simple case G l MR = MMR/ S \ ars , where / s \ ars is the (constant after t- m ) mass fraction of the 
it/i-element in evolved (mass-losing) stars. The SNIa enrichment term G SNIa = iVsNia2/sNia> 
where 2/g NIa is the ith element SNIa mass yield. The SNIa rate may be characterized as 



iV S NT a ~ BL R (t/t n(mr ) & with B = 1.6 10 13 fl sNTa y r 1 L^- # S Nia = 1 corresponds to the 



rate of ICappe 



the models of 



laro et al . f 199 71). ICiotti et all (119911) adopt 



7 SNIa 

1.1 and # S Ni a 



1, while in 



Ruiter. Belczynski. Sz Fryer! ( 120091 ) substantially lower rates, and a somewhat 
steeper decline, are derived from binary population synthesis. In some models we assume 
that there is an outflow that is driven by SNIa and given by M out = i^out-^SNia-^iSM, where 
If out is a dimensionless constant. 



A. 2. Special Cases with Exact Solutions 

It proves convenient to introduce the following dimensionless variables: r = t/t now , fi 

^ISMA^ISM^now); 0* = /iSM //stars) ^ = Mt now / Mi SM (t now ) , Y = G*t now /(M ISM (t now )/* tars 
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Equations (Al) and (A2) become 

dfliSM 

dr 

and 



= ^MR - VWt, (A3) 



^(^ism^ism) i i n , , x, ,i / A/| x 
J t = 7sNIa + 7MR - I 1 + Ooutj^out^lSM, ( A4 ) 

with = 1 and the following boundary conditions at r in = t- m /t now : /i in = MisM^inV^iSM^now) 
and 0i SM ,in = /ism (*in) //stars- The source terms are ^ M r = i^ a , 7^ R = ^mr, 7sNia = 
e i ?] r- /3 , with^ = A(M 1SM (t now ) / L B )t now , r] Q = B(M ISM (t now )/L B )t now , and e l = y l sm J /* tars ; 
for SNIa-driven galactic winds tp out = Kr/oT" 13 /i, where k = K out Mi SM (t now ) . Q % = r] e % /ip = 
(B/A)e l and W = k,t] /^ = (B/A)k are convenient dimensionless quantities that factor into 
determining the importance of SNIa metal enrichment and outflows, respectively. 

There are two classes of simple analytic solutions to equations (A3) and (A4) - closed 
box models where M ou t = 0, and "tracking" solutions where a = (3. Subsets of the latter 
include steady-state solutions (only exactly possible if a — j3), and solutions where the ISM 
mass is in a steady state, but its metallicity is not (i.e., abundances are evolving much more 
quickly than the total mass). 

For closed box models, the dimensionless abundances at the present day are given by 

<«1) = 1 " Ml - 0isM,in) + (TT^a - rL- ), (A5) 
where ji- in or r in may be eliminated using 

^ = 1-7^(1-^). (A6) 
{I -a) 

This takes a particularly simple form for the reasonable initial conditions 4>\ SM in = 1 (iden- 
tical ISM and stellar metallicities) and/or /i in = (galaxy evacuated of gas at the close of 
the star formation era). 

A steady state is possible if a — j3. The dimensionless ISM mass fraction at all times 
in this case is given by 

(1+g*) 

0ISM = n V ( A7 ) 

(1 + Oout) 
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We also consider the general a — /3 case with SNIa-driven outflow and b out = 0. In this 
case, the solution is 

0ism(1) = > (A8) 

where W and fi in are related to r in via 

/i in = W- 1 {1 - (1 - W)exp[W^ {\ - rl- a )/{\ - a)]} (A9) 



A. 3. General Features and Application 

All of the simple models above have solutions of the form 

ti S M(l) = r i + sqj,i.e- (A10) 

/bm(!) = r Vstars + ^SNIa, ( A11 ) 

where q = B/A is proportional to the present-day ratio of SNIa to mass return rate. For 
the simplest steady state models (r l — s — 1) mass return enriches the ISM to the stellar 
metallicity (first term), and SNIa further enhances the abundances according to the ratio of 
supernova to mass-loss rates and the SNIa yield. In the other models considered above r 
represents dilution or amplification of stellar mass return due to initial conditions for closed 
box models, biased outflow in the steady-state b out ^ case, and the interplay between 
initial conditions and outflow for the tracking solutions. The factor sq represents the SNIa 
enrichment and depends on the ratio of supernova to mass-loss rates with these same causes of 
dilution/amplification for the steady-state 6 out ^ and tracking solutions, and modification 
due to the mismatch in the SNIa and stellar mass loss time dependencies for the close box 
case. 

In particular (see Table 3) r l — (1 + & ou t) _1 for steady state, r % = 1 — /j, in (l — (p\ SM in ) for 
closed box, and r % = [(1 — fi in ) + A*in0isM in(l — ~ WMin) f° r tracking solutions. The 

SNIa parameter s = (1 + 6 ou t) _1 for steady state, s = ip (l — _/3 )/(l — (3) for closed box, 
and s — (1 — /i in )/(l — Wfi in ) for tracking solutions. 

One can relate the ISM abundances of the ith and jth elements, eliminating the factor 
s, via 

/iSM(l) = r /stars + (/iSm(1) — r /stars) (A12) 

2/SNIa 
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or 



/] 



ISM^ 



fi fj 
J stars J si 



stars 



1SM\ 



fj fJ 
J stars J If- 



■r + I 1 — r 



ISM \ 




(A13) 



where we have made the simplifying assumption that r l = r is the same for all elements (al- 
ready implied for steady state solutions by adopting constant 6 out , and otherwise equivalent 
to constant <j)\ SM in - i.e., an element-independent initial ratio of ISM-to-stellar abundances). 

In the approach we adopt here, we calculate the ISM abundance pattern for a given 
stellar abundance pattern (corresponding to a particular [a/Fe] ratio; see below) and ISM- 
to-stellar Fe abundance, for a range of values of r. One then compares the pattern with the 
NGC 4472 observations to see if consistent scenarios emerge and to identify any anomalies. 
The parameter r may then be interpreted post facto under the different classes of models, 
as can s for values of the stellar and ISM Fe abundances (and not just the ratio s ). W e 
consider the SNIa and SNII+HN yields from lNomoto et al.l ( 119971 ); iKobayashi et al.l ( 120061 ). 
and use these to set the stellar abundance pattern for a given value of [a/Fe] s tars, where 
[a/Fe] is defined as follows. We calculate the ratio, with respect to the Fe mass fraction, 
of the sum of the mass fractions of all elements under cons i derat ion lighter than Fe, divide 
by the same ratio calculated using the iGrevesse Sauvall (119981 ) solar standard, and take 
the logarithm. For a given set of yields, [a/Fe] corresponds to a particular number ratio of 



SNIa to SNII that contributed to establishing the stellar abundances. The IKobayashi et al. 



( 120061 ) SNII+HN yields are averaged over a Salpeter IMF with maximum mass 50 M , and 



the hypernova fraction set to match the abundances of low metallicity stars in the Milky 
Way halo. Because of the rapi d buildup of metals in e lliptical galaxies, we adopt the solar 
metallicity SNII+HN yields of IKobayashi et al.l (120061 ) . Selection of a different metallicity 
would have little effect on the final stellar abundance pattern and corresponding value of 
[a/Fe] (compared to the uncertainty in the yields themselves; Gibson et al. 1997), and 
primarily alter the corresponding relative numbers of SNII and SNIa. The effects of different 
yield sets and IMFs is beyond the scope of the present work, but something we will consider 
in future investigations. The same set of yields are used for SNIa that directly enrich the 
ISM, and those that enrich the stars. 



Table 3. Definitions for r and s in Different Models 



Model 


r i 


S 


Steady state 


(l + boutT 1 


(l+fcoutT 1 


Closed Box 


i - Wn (i - 4>\su,iJ 




Tracking 




(1 - Wn )/(1 - Witta) 
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As it turns out, numerical experiments demonstrate that s ~ r for most reasonable 
models. This is exactly true for the steady state models. For closed box models with 
^iSM,in = °) ^ deviates by < 20% for 6 MR = 1, Mi SM (t now ) / L B ~ 0.1, and < < 2. The 
dependence on (p\ SM in is flat until one gets to the regime where the solution corresponds to 
the ISM Fe abundance being substellar until very recently and then suddenly being enriched 
to the observed superstellar value by a concentrated period of SNIa enrichment. Similarly, 
for tracking solutions s differs significantly from r only in cases where, at a recent epoch 
the galaxy was filled with a large mass of low-abundance gas. In this scenario, the galaxy 
must be experiencing a powerful outflow that has ejected most of the low metallicity gas 
and replenished it with SNIa-enhanced stellar mass return. For example, using our standard 
parameters, to obtain s = 0.5 for r = 1, one requires the galaxy have 6.4 times its present 
ISM mass at r in = 0.7, and with < 0.08 solar Fe abundance. This implies that, for relatively 
quiescent elliptical galaxies, for a given value of the dilution factor r the offset between the 
stellar and ISM Fe abundances (and the SNIa Fe yield) can be robustly used to estimate the 
SNIa to mass return rate. If the SNIa and mass loss rates have different time- dependencies 
the ratio will correspond to a ratio of time- averaged rates, but for standard parameters this 
does not drastically differ from the ratio at the present day. 
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